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We derive the generalized Boltzmann equation (GBE) near equilibrium from the Kadanoff-Baym 
equation for quark excitation with ultrasoft momentum (~ g 2 T, g: coupling constant, T: tem- 
perature) in quantum chromodynamics (QCD) at extremely high T, and show that the equation 
is equivalent to the self-consistent equation derived in the resummed perturbation scheme used to 
analyze the quark propagator. We obtain the expressions of the dispersion relation, the damping 
rate, and the strength of a quark excitation with ultrasoft momentum by solving the GBE. We also 
show that the GBE enables us to obtain the equation determining the n-point function containing 
a pair of quarks and (n — 2) gluon external lines whose momenta are ultrasoft. 
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I. INTRODUCTION 



Theories containing fermion and boson such as the 
Yukawa model, quantum electrodynamics (QED), and es- 
pecially quantum chromodynamics (QCD) suggest that 
the plasmas described by these theories including QED 
plasma and the quark-gluon plasma [1] at so high tem- 
perature (T) that the particle masses are negligible are 
multi-scale systems. For example, the average inter- 
particle distance is of order T _1 , the Dcbyc screening 
length, (gT)~ x (g: coupling constant) [2], and the mean 
free path of the quark and the gluon, (g 2 T) _1 [3, 4], re- 
spectively. These momentum scales, T, gT, and g 2 T, 
are called hard, soft, ultrasoft, respectively. Due to the 
existence of these multiple scales, analyzing the spec- 
tral properties of collective excitations at finite temper- 
ature in these fermion-boson systems is a quite nontriv- 
ial task even at weak coupling (g <C 1) regime in gen- 
eral. In the soft scale, the well-established perturbation 
theory known as the hard thermal loop (HTL) approxi- 
mation [5] is applicable, and the analysis using that ap- 
proximation suggests the existence of the plasmon [2], 
which is a bosonic excitation and well-known in nonrcla- 
tivistic plasma. In addition to the bosonic excitation, a 
fermionic excitation called plasmino [6] is also suggested 
to exist, owing to the fact that the particle masses are 
negligible compared with T. By contrast, the HTL ap- 
proximation can not be used in the ultrasoft scale due 
to an infrared singularity known as pinch singularity, 
and thus a novel resummed perturbation theory, which 
regularizes the singularity, needs to be developed to ob- 
tain the correct leading-order result [4, 7, 8]. Due to 
the same reason, the resummed perturbation theory is 
also necessary in the calculation of the bosonic propaga- 
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tor [9] and the transport coefficient [10-13]. The anal- 
ysis using the resummed perturbation revealed the exis- 
tence of a novel quark excitation (ultrasoft fermion mode) 
whose momentum is ultrasoft [4] in addition to the plas- 
mino [6] , and the expression of the pole position and the 
strength of that excitation was obtained by solving the 
self-consistent equation derived in the resummed pertur- 
bation, in the Yukawa model and QED [7]. Wc note 
that results which suggest the existence of the ultrasoft 
fermion mode were also obtained by other analyses, in 
which the Schwinger-Dyson equation [14], the Nambu- 
Jona-Lasino model [15], and effective models in which 
the boson has finite bare mass [16] are used. On the 
other hand, in QCD, the self-consistent equation has not 
been solved, and thus the expressions of the pole position 
and the strength of the ultrasoft fermion mode have not 
been obtained. In this paper, we extend the analysis in 
the Yukawa model and QED [7] to QCD, and obtain the 
expression of the pole position and the strength of that 
excitation. 

Now we introduce another point of view: each per- 
turbation scheme is equivalent to different kinetic equa- 
tion near equilibrium. In fact, the HTL analysis of the 
soft boson (fermion) propagator is equivalent to the usual 
(generalized) collisionless kinetic equation called (gener- 
alized) Vlasov equation [17]. Here the generalized Vlasov 
equation means the collisionless kinetic equation in which 
there is a fermionic background field in place of bosonic 
one, in contrast to the usual Vlasov equation case. On the 
other hand, it was shown that the resummation scheme 
used in the analysis of the gluon propagator whose mo- 
mentum is ultrasoft, is equivalent to the Boltzmann equa- 
tion [9]. It reflects the fact that the interaction among 
the particles can not be neglected when the space-time 
scale is comparable with or much larger than the mean 
free path, which is of order (g 2 T) _1 in the case of the 
hard quark and gluon [3, 4]. We note that this discus- 
sion explains the reason why the resummed perturbation, 
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which takes into account the interaction effect among the 
hard particles, is necessary in the analysis of the ultra- 
soft momentum region. The resummaton scheme used 
in the analysis of the ultrasoft fermion propagator [4, 7] 
is also interpreted as a generalized version of the Boltz- 
mann equation (GBE) near equilibrium in the case of the 
Yukawa model and QED [18]. In this paper, we derive 
the GBE that is equivalent to the self-consistent equation 
derived in the resummed perturbation [4] in QCD. 

The aim of this paper is twofold: One is to obtain the 
expression of the dispersion relation, the damping rate, 
and the strength of the ultrasoft fermion mode in QCD. 
The other is, to derive the GBE that is equivalent to the 
self-consistent equation obtained in the resummed per- 
turbation [4] from the Kadanoff-baym equation [19, 20] 
in QCD. The derivation helps us to establish the foun- 
dation of the resummed perturbation scheme [4], and 
to obtain the kinetic interpretation of the procedure of 
that scheme. As a by-product, the derivation using the 
Kadanoff-Baym equation enables us to evaluate also the 
nonlinear response that is caused by the average gluon 
field, not only the linear response caused by the quark 
average field [9, 17, 19]. Due to this merit, we can ana- 
lyze the n-point function whose external lines are (n — 2) 
ultrasoft gluons and a pair of ultrasoft quark, not only 
the ultrasoft quark propagator. This paper is an general- 
ization of Ref. [7, 18], from the Yukawa model and QED 
to QCD, and also an extension of Ref. [9], from the gluon 
propagator to the quark propagator. Since there is self- 
coupling of the gluon and the color structure of the off- 
diagonal propagator, which is an important quantity in 
our analysis and will be introduced later, is complicated 
in general in QCD, the extension to QCD is nontrivial 
and thus it deserves a paper. 

This paper is organized as follows: We derive the GBE 
near equilibrium that is equivalent to the self-consistent 
equation obtained in the resummed perturbation the- 
ory [4, 7] in linear response regime in Sec. II, in a similar 
method to the case of the Yukawa model and QED [18]. 
We also show that the Ward-Takahashi identity is sat- 
isfied, and that the identity can be derived from the 
conservation of the color current. In Sec. Ill, we ana- 
lyze the properties of the ultrasoft fermion mode such as 
dispersion relation, damping rate, and residue in QCD, 
by solving the GBE in some momentum region as in the 
case of the Yukawa model and QED [7] . Section IV is de- 
voted to obtaining the equation determining the n-point 
function whose external lines are a pair of quarks and 
(n — 2) gluons with ultrasoft momenta. We summarize 
our paper and give concluding remarks in Sec. V. In Ap- 
pendix A, we show that the longitudinal component of 
the off-diagonal propagator introduced in Sec. II is much 
smaller than the transverse component of that quantity. 
We give a detailed derivation of Eqs. (2.51) and (2.64) 
in Appendix B. Appendix C is devoted to derivation of 
Eq. (3.1). 



II. GENERALIZED BOLTZMANN EQUATION 

In this section, we derive the GBE that is equiva- 
lent to the self-consistent equation obtained in the re- 
summed perturbation theory [4], from the Kadanoff- 
Baym equation [19, 20] by using the background field 
gauge method [21], which enables us to obtain the resul- 
tant equations in a gauge-covariant way. We also check 
that the resultant kinetic equation satisfies the Ward- 
Takahashi identity, which is a necessary condition of the 
SU(N) gauge invariance. 



A. Background field gauge method 

Let us introduce the background field gauge 
method [19, 21]. In this method, the following gener- 
ating functional is used: 

Z\j, T), fj; A, V,W} = J [Va] [V$] [2ty] [D(] [DC] 

xexp(i f d 4 x(C[A*+a a ^,V + ij,W + ^}+£ F p (2.1) 
J c 

-0> OM +^ + ^))). 

where the time integral is defined using the closed-time- 
path formalism [19, 22]: J c d 4 x = J c dx° J d 3 x, where 
J c dx° is the integral along the contour C, which is de- 
fined in Fig. 2 and consists of C + , C~ , and C°. Here tf 
is larger than any time we consider. A° and are the 

vector fields, & (&) and ip (ip) the (anti-) spinor fields, ( 
(C) the (anti-) ghost field, the external current, and 77 
(rj) the external (anti-) fermion source, respectively. 
The Lagrangian of QCD is defined as 

Nf _ 

C[a,i>,i>] =iJ2^Pl<Ai>j - /"^HCW. (2-2) 
3=1 

where F£ v = d IJ .a1—d v a'^—gf ahc a h ll a^ / is the field strength, 
D^[a] = d^+iga^t 11 the covariant derivative in the funda- 
mental representation, t a (a = 1, ...TV 2 — 1) the generator 
of SU (N) group in the fundamental representation, f abc 
the structure constant of SU(N) group, respectively. We 
note that Nf, the flavor number, and N, the color num- 
ber, are not specified in this paper, because we want to 
see how the expressions of the properties of the ultrasoft 
mode, which will be analyzed in Sec. Ill, depend on N 
and Nf. In the real world, N = 3. Since we are consid- 
ering the case that the quark current mass is negligible 
in every flavor, we drop the index for the flavor from now 
on. 

We adopt the temporal gauge fixing, which was used 
in the analysis based on the resummed perturbation [4]. 
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FIG. 1: The correspondence between the diagrammatic methods and the kinetic equations, and the fermionic modes that result 
from the analysis using the two methods. The vertical axis denote the energy scale. 
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FIG. 2: The contour C in the complex x° plane. 

In this gauge-fixing, the Faddeev-Popov term is 
Cfp = _ A _ f ( * S* - r d \A + a){) C b , 



(2.3) 



where the gauge- fixing function is G a [a] = ag with A 
oo, which is equivalent to the following constraint: 



a a = 0. 



(2.4) 



We note that the possible gauge-fixing dependence needs 
to be checked apart from the covariance with respect to 
the background gauge transformation, which will be in- 
troduced later. 

In the background field gauge method, we impose the 
following conditions: 



K) = W = <V> = o. 



(2.5) 



A^, \P become equal to the gluon and the (anti-) 
quark average field after imposing these conditions [21]. 
We note that the action in Eq. (2.1) is invariant under 



the following transformation [19]: 

ip(x) — > h(x)ip{x), ?p(x) —> ip(x)h^(x), 
A a ^(x)t a -> h{x)Al(x)t a h\x) - l -hd^h\x), 

a%{x)i? -»• h(x)a«(x)t a h\x), (2 6) 

j;(x)t a -> h(x)%{x)rht{x), 

r](x) — > h(x)rj(x), rj(x) — > rj{x)h/{x), 
C(x)t a -> h{x)C{x)t a h\x), 
T{x)t a ^h\x)C{x)t a h{x), 

where h(x) = exp[i8 a (x)t a ]. 



B. Derivation 

We consider the following situation: The system is 
at equilibrium whose temperature is T before the initial 
time to. Then, external quark, anti-quark, and gluon ex- 
ternal sources disturb the system, and hence the system 
becomes non-equilibrium state. In this paper, we focus 
on the case where the external fermion source is so weak 
that we need to retain only the contributions that are in 
the linear order in the magnitude of the fermionic aver- 
age field \P to the induced fermion source, which will be 
introduced later. Due to the linear response theory, anal- 
ysis in such situation is equivalent to the computation 
of the fermion propagator with ultrasoft momentum at 
equilibrium. We note however, that the induced fermion 
source contains the all order contributions in A°; in our 
approximation, in which we have a good machinery to 
compute the n-point function whose external lines are 
two quarks and (n — 2) gluons and whose all momenta 
are ultrasoft, not only the ultrasoft fermion propagator. 
In this sense, we are going to analyze the region which is 
beyond the linear response regime. The derivation will 
be performed in a similar way to that in Ref. [9, 18]: We 
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will apply the gradient expansion and the weak coupling 
approximation to the Kadanoff-Baym equation [19, 20]. 
We note that both of them are justified by the smallncss 
of the coupling constant, as we will see later. 

By performing a infinitesimal variation of the integral 
variable in Eq. (2.1), we get the following equations: 



{ip x [A + a](& + ip)(x)) =r 1 (x). 



(2.7) 



(2.11) 



((D U [A + a]) ab F\[A + a\- g{(M> + <P)t a ^(lP + V>)> 

~f abc 9o»(C b ( c )=J a ,(x). (2.8) 

Here [a] = + iga^T a is the covariant derivative in 
the adjoint representation, where (T a )i, c = —if abc is the 
generator of SU (N) group in the adjoint representation. 
By imposing Eq. (2.5), we get the following equations of 
motion for the average fields: 

ip x [Ap(x) = T)(x) + 7] ind (x), 

(2.9) 

D? x [A]F b ^[A](x) - gW(x)t a j^(x) = f^x) + jZ(x). 

(2.10) 

Here ?7i n d and j°^ A are the induced fermion source and 
the induced color current defined as 

?7ind(^) = gt a (4° '(x)ip(x)) , 
jZ(x) = g^t a r^ + f abc go»(C b C) 

+ gf abc r^(ald xa c) + g 2 nZ X Al(a c p a d x ) 
+ g 2 f abc f cde (a b l/ a dp a e "), (2.12) 

where we have introduced r^ pXv = 2g pp g Xu - g pX g pv ~ 
g pv g pX and 

KZ X = \[f eab r d (g^9 uX - g^g vp ) 

_|_ jeac jedb^X gpv _ giiu gXp-j (2-13) 

+ f ead f ebc {g fJ ' l/ g pX — 9 pp 9 Xu )]- 

It was shown [19] that Eqs. (2.9) and (2.10) transform 
under the gauge transformation given in Eq. (2.6) in a 
covariant way. 

By differentiating Eq. (2.7) with respect to j ap {y) and 
Eq. (2.8) with respect to rj(y), we get the following equa- 
tions: 



ip x [A]K a Jx, y) - g>fD?Jx, y)t b *(x) = * 



. Sr] ind {x) 
Sj pa (y) ' 



(2.14) 



{{D 2 [A]) ab g^ - {b p [A)b v [A)) ab + 2igF cpv [A\T c ah ) x 



xK(y,x) 



(2.15) 



g (#(x)t a rWxMy)) + S(y,x)t a ^W(x)) 



MM 



where we have introduced the following propagators: 

Df v (x,y) = (T c a»(x)al(y)), (2.16) 
S(x,y) = (Tc1>{x$(y)), (2.17) 
K^x,y) = (T c ^{x)a^y)). (2.18) 

Here Tc means the path-ordered product whose path is 
C: 

D^( X ,y) = 6c(x ,y°)D>f(x,y) 



+ e c (y a : x )D<? b (x,y), 



(2.19) 

S(x, y) = 9 c (x°, y°)S> (x, y) - 6 c (y°, x°)S<(x, y), 

(2.20) 

K;(x, y) = c (x°, y°)K> a (x, y) + 6 c (y°, x°)K< a (x, y). 

(2.21) 

where we have introduced the step-function along the 
contour C, 9c(x,y), and the following functions: 



D>? b (x,y)^ 


(o-(x)^(y)), 


(2.22) 


D^ ab (x,y)^ 


(aUvK(x)), 


(2.23) 


S>(x,y) = 


(ip(x)ip(y)), 


(2.24) 


S < (x,y) = 


(ip{y)tp(x)), 


(2.25) 


Kl a {x,y)^ 


mxK(y)), 


(2.26) 


K< a (x,y) = 




(2.27) 



In fact, K> a (x,y) coincides with K< a (x,y) in our ap- 
proximation as in the case of the Yukawa model and 
QED [18], so we simply write these two functions as 
K°i(x,y) from now on. We call K°i 11 off-diagonal prop- 
agator" since that quantity is the propagator between 
the quark and the gluon [18]. We note that K°i vanishes 
at equilibrium. The off-diagonal propagator is the most 
important propagator in our analysis, which can be seen 
from the discussion in Sec. II D. The transformation form 
of these propagators under the gauge transformation is 
given in Ref. [19]. 

From Eq. (2.4), we get 



D% b (x,y)=D<$ ) (x,y) = 0, 
Kg(x,y)=0. 



(2.28) 
(2.29) 



Because of these constraints, we do not have to deal with 
the temporal components of these propagators. 

By setting x° £ C + and y° £ C~ , and interchanging x 
and y in Eq. (2.15), we get 

iP x [A]K?(x,y) - g 7 Wf ba (x,y)t b *(x) = 

(2.30) 

((D 2 {A]r b g v - (D i [A]D j [A]) ab + 2igF^[A]T c ab ) y 
xK b (x,y) 



Sr)(y) 



gS<(x,y)t a fV(y) = t- 



Srj(x) 



(2.31) 
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where we have neglected 'P(y)(ip(y)ip(x)) because 
(ip(y)'ip(x)) contains more than one average quark field. 
The induced source terms arc evaluated as follows [18]: 

Sf a (y) 

= f d i z{E{x ) z)K1{z 1 y) + S h ^x 1 z)D)t{z,y)) 
Jc 

POO 

i I d 4 z 

(2.32) 



E m \x,z)Df ba {z,y)), 



5r)(x) 

= f d 4 z(n aM i(z,y)K»(x,z)+~ m (z,y)S(x,z)) 
Jc 

/oo 
d 4 z(n Aab ^(z,y)K b (x,z) 
-oo 

-S Rai (z,y)S<(x,z)). 



(2.33) 



Here £ is the quark self-energy, Pi^ b the gluon self- 
energy, and the off-diagonal self-energy [18], respec- 
tively, and they are decomposed as follows: 

S(x, y) = 6 c (x°, y°)S > (x, y) - 9 c (y°, x°)S<(x, y), 

(2.34) 

nf v {x,y)=ec{x\y°)Ilf>{x,y) 

+ d c (y ,x )n?<(x,y), (2.35) 

Sftx, y) = 6 c {x\ y°)~«> (x, y) + 6 c {y°, x°)^< (x, y). 

(2.36) 

We have also introduced the retarded quark self-energy 
S R , the advanced gluon self-energy n A " b , and the re- 
tarded off-diagonal self-energy S Ra , which are defined as 

E R (x, y) = i9(x°, y^l^ix, y) + E<{x, y)}, (2.37) 

n Aab (x,y) = -i6(y ,x°)[n; b >(x, y ) - nf<(x, y )\, 

(2.38) 

S Ra (x,y) = ie(x°,y°)[S;>(x,y) - S a <{x,y)}. (2.39) 

From Eqs. (2.30)~(2.33), we have 

ip x [A]K°(x, y) - gi Wf ba (x, y)t b *(x) (2.40) 

#z(e r (x, z)K?(z, y) + S m i (x, z)D< ba (z, y)) , 

((D 2 [A]) ab g i3 - (D' l [A]D j [A]) ab + 2igF cl] [A]T^ b ) y 
xK b (x lV ) 

/OO 
d ^n Aa ^(z,y)K b (x,z) 
-oo 

-5 Rai (z,y)S<(x,z)). (2.41) 



Now we perform the Wigner transformation [19] , which 
is defined as 



f(k,X) = J d 



(2.42) 



where s = x— y, X = (x+y) /2, and f(x, y) is an arbitrary 
function. After this manipulation, Eqs. (2.40) and (2.41) 
become 



(X)t b K?(k,X) 



91 W< ba (k,X)t b <P(X) 



= £ H (k,X)K?(k,X) + S Hb \k,X)D< ba (k,X) 
{[(-k 2 + ik ■ d x )S ab - 2gk ■ A c (X)(T c ) ab ]g k i 



(2.43) 



-k k k j + -(d x k j +k k d J x ) 



tub 



-{ 

- g{T c ) ab (k k A CJ {X) + A ck (X)k j )^K b (k,X) 
+ gS < (k,X)t a j k <P{X) = II Aabk 3{k,X)K b (k,X) 

- ~ Rak (k,X)S < (k,X). 



(2.44) 



Here we have used the gradient expansion [9, 17] by as- 
suming k ;§> dx, as in Ref. [18], and set i — > k in the latter 
equation. This assumption is justified [18] by using the 
fact that k ~ T and dx ~ g 2 T, which is valid since we 
are focusing on the ultrasoft energy region. By multiply- 
ing Eq. (2.43) by (ft + i$x/2 - gA b (X)t b - Z R (k,X)), 
multiplying Eq. (2.44) by the projection operator into 
the transverse component P^ k (k), which is defined as 

= 9tii9v](.&ij ~ hkj) with ki = fcj/|k|, and sub- 
tracting the latter from the former, we obtain 

(2ik ■ d x - 2gk ■ A b (X)t b - {ft, E R (k, X)}) K?(k, X) 
-2gk-A c (T c ) ab K b (k,X) 



+ pT k (k)n Aabk Hk,x)K b (k,x) 

= g{m ba (k 1 X) + Pj l {k)5 ab S < {k))r b i{k,X), 



(2.45) 



with gt^(k, X) = gt a j^(X) + S Ra (k, X). Here we have 
neglected the longitudinal component of Kf(k, X), which 
is much smaller than the transverse one as is shown in 
Appendix A, and used Pf k (k)K ak (k,X) ~ —Kf(k, X). 

The diagonal propagators in the right-hand side can 
be replaced by that in the free limit at equilibrium since 
we are considering the case where the system is near the 
equilibrium so that the induced fcrmion source contains 
only one average fermion field (>P) [18]: 



S a< (k) = ftp Q {k)n F (k°), 
D°< ab {k) = 6 ab p°(k)n B (k a )pT(k), 



(2.46) 
(2.47) 



where p°(k) = 27rsgn(fc°)<5(fc 2 ) is the free spectral func- 
tion and n F (k°) = (e^ + l)- 1 (n B (k°) = (e k °/ T -l)- 1 ) 
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is the fermion (boson) distribution function at equilib- 
rium. We note that only k satisfying k 2 = contributes 
to the fermion induced source since 5(k 2 ) appears in the 
right-hand side of Eq. (2.45). 

Now let us see that the self-energy terms in Eq. (2.45) 
is not negligible in comparison with the first term in the 
left-hand side, i. e., 2ik ■ dxKf (fc, X). First, we evaluate 
the diagonal self-energy terms. The dominant parts of 
the diagonal self-energies are those at the leading order 
and at equilibrium, whose diagrams are shown in Figs. 3 
and 4 and whose expressions arc given as follows: 



{%E R ^\k)} =m 2 g -2i( g k°, 
P^i(k)n A ^ abau (k) = -5 ab P? v (k)(m, 



(2.48) 



2i( g k°), 



(2.49) 



where m q = gT^/Cf/2 (m g = gT^/(N + N f /2)/6) is the 
quark (gluon) asymptotic thermal mass [23], with Cf = 
(N 2 ~1)/(2N). ( q ((g) is the quark (gluon) damping rate, 
which is of order g 2 T\n(l/ g) [3]. We have used the fact 
that k ~ T and k 2 = 0. Since Eqs. (2.48) and (2.49) are 
of order g 2 T 2 , we confirm that the diagonal self-energy 
terms in the left-hand side of Eq. (2.45) have the same 
order of magnitude as that of 2ik ■ dxKf(k, X) in the 
left-hand side because k ~ T and dx ~ g 2 T. For this 
reason, we can not neglect the diagonal self-energies in 
that equation, in contrast to the soft case (dx ~ gT) [17]. 

As in the diagonal self-energies case, the off-diagonal 
self-energy need to be taken into account. The expression 
of the off-diagonal self-energy at the leading order, whose 
diagrams are shown in Fig. 5, is given by the following 
equation: 



E: Ha »(k,x) = 
(M 



■jH b S 0R (k + k')^t a K b Ak', X) 



- g2 J ^ k ^T d (k'-k) 

x if abc [g lj (2k' - kf + gM(-k - k') 1 
+ g^(2k-k'y]K b (k',X), 



(2.50) 



- k 2 ] and D° l Rcd (k) 



where S 0R (k) = -$/[(k + 
-S cd {P^(k)/[(k +ie) 2 -k 2 }+kikk/(kQ+ie) 2 }am the free 
quark and gluon retarded propagator at equilibrium, re- 
spectively. By making an order estimate of Eq. (2.50), we 
see that the off-diagonal self-energy term in Eq. (2.45) has 
the same order of magnitude as that of 2ik ■ dxKf(k, X), 
and thus that term can not be neglected. 

By inserting Eqs. (2.46)~(2.50) into Eq. (2.45), we ob- 
tain the following equation (see Appendix B for a detailed 



derivation) 



[2ik ■ d x - 2gk ■ A b (X)t b + 2i(k° + 5m 2 ) Kf(k, X) 
-2gk-A c (T c ) ab K b (k,X) 
= k : p Q (k)pT(k)(n B (k°)+n F (k )) 



(2nY 



k ■ k' 



(2.51) 



where we have introduced 



r 2 2 

dm = m„ — m 



g *T 2 + (2.52) 

y 1 24 8N 12 ' 



and C = Cq + C g , and used the relation P^.(k)K ak (k, X) ~ 
—Kf(k,X) again. We note that this equation trans- 
forms in a covariant way under the gauge transformation 
Eq. (2.6), which can be confirmed by using the transfor- 
mation property of Kf(k, X) [19], 



K?(k,X) -> h(X)K?(k,X)h ab (X) 



(2.53) 



where h(x) = exp[i9 a (x)T a }. The diagram that corre- 
sponds to Eq. (2.51) is shown in Fig. 6. From Eq. (2.51), 
we get the following equation for Ki(k, X) = t a Kf(k, X): 

(2ik ■ d x - 2gk ■ A a (X)t a + 2i(k° + 5m 2 ) Ki(k,X) 

= gC f ^p°(k)P^(k)(n B (k ) + n F (k°)) 

( r d 4 k' k k ^y j + ^ k k' j 

X h^(X)+g / ° * k \ + y K k (k',X] 



(2tt)* 



k ■ k' 



(2.54) 



From this equation, we confirm that Ki(k, X) transforms 
under the gauge transformation Eq. (2.6) in the same way 
as &(X) in contrast to Kf{k,X). 

By introducing A^±(k, X), which is de- 
fined by K^(k,X) ee 2Tr5(k 2 )[6(k )A fl+ (k,X) + 
6(-k°)A^(-k,X)}, we arrive at the following GBE 1 : 



2iv ■ D x [A] ± 



5m 



2<U±(k,X) 



2gC fl 6[n B (\k\) + n F (\k\)}#(X) 

-g 2 C fl ^[n F (\k\)+n B (\k\)]P^(v) 

v V- / d 3 k> 1 , r ,, 

X JLJ ( 27r )3 2 |k'| \k\\k>\w> A ^ k > X >> 

(2.55) 



1 In the case of the analysis of the lightlike momentum instead 
of the ultrasoft one, a equation which is similar to this equation 
was obtained before [25] . However, the equation in Rcf. [25] does 
not contain all the leading contributions: The equation does not 
contain the third term in the left-hand side and the second term 
in the right-hand side in Eq. (2.55). 




FIG. 3: The quark retarded self-energy S R ^ e(1 \k) whose mo- 
mentum is hard in the leading order. The solid line is the 
quark propagator and the cirly line is the gluon propagator. 
We note that the gluon propagator in this figure is the HTL- 
resummed one [24] , which results the anomalously large imag- 
inary part (C„) of r fl(cq) (fc). 
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FIG. 5: The off-diagonal self-energy S^(k,X) whose momen- 
tum is hard in the leading order. The propagator that is 
composed of the solid line and the curly line with the black 
blob is the off-diagonal propagator. The other notations are 
the same as Fig. 3. 
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FIG. 4: The gluon advanced self-energy n A ^ ab ^(k) whose 
momentum is hard in the leading order. The notations are 
the same as Fig. 3. The gluon propagators in the second 
diagram are the HTL-resummed ones [24], which result the 
anomalously large imaginary part (C s ) of i7 A(cq)ai> '" / (fc). The 
diagram which has ghost propagators is omitted since its con- 
tribution to the transverse sector of 77 A(eq)a6A "'(fc) is zero. 



where v>* = (l,k) and v ,fi = (l,k'), with k = k/|k|. 
Since the structure of this equation is the same as that 
in QED [18], the interpretation of each term in that equa- 
tion is the same as in the QED case [18]. 



C. Linear response regime 

When A£ = 0, the color structure of the off-diagonal 
propagator becomes simple. In this case, Eq. (2.51) be- 
comes 



(2ik ■ d x + 2iC,k° + 5m 2 ) K?(k, X) 



o°(k)pT(k)(n B (k°)+n F (k )) 



(2tt) 4 k ■ k' 



(2.56) 



This equation tells us that Kf has the following color 
structure: 



K?(k,X)=t 



a K i(k,X) 



(2.57) 



In the present case, Ki(k, X) has the same color structure 
as that of \P(X). For later use, we point out that has 



also the same color structure as Kf . By using Eq. (2.57), 
Eq. (2.56) becomes 

(2ik ■ d x + 2i(k° + 5m 2 ) K^k, X) 



gC f ^ P °(k)P^(k)(n B (k ) + n F (k )) 
d 4 k' k k -y j +j k k' j 



x v>*(X) + g 



(2nY 



k-k' 



K k (k',X] 



(2.58) 



D. Equivalence between generalized Boltzmann 
equation and resummed perturbation 

Let us show that the GBE in A^(X) = case is equiv- 
alent to the self-consistent equation derived in the re- 
summed perturbation [4] . For this purpose, we write the 
off-diagonal propagator in terms of the off-diagonal self- 
energy using Eq. (2.58): 



(2k ■ p + 2i(k° + 5m 2 ) K^k, p) 

\k)P^(k)(n B (k°) + n F (k°))fj(k,p). 



(2.59) 



Here we have performed the Fourier transformation as 
f(p) ee / d 4 X/(2n) 4 e i P- x f(X), with f(X) is an arbitrary 
function, and introduced r j (k,p), which is defined as 
r a i(k,p) = t a fi{k,p)/C f . Since Sf has the same color 
structure as that of so does pi . Now we can obtain 
the self-consistent equation in terms of the off-diagonal 
self-energy, by inserting Eq. (2.59) into Eq. (B3): 



r l (k,p) 



d 4 k' wy + ^k' 1 



C f ' ' * J (2tt) 4 k-k' 
Vp°(k')Pj k {k'){n B {k' G ) + n F (k' )) 
2k' ■ p + 2i(k'° + 8m 2 



r k (k', P ), 



(2.60) 



where we have introduced Pi(k,p) by Fi(k,p) = 
Fi(k,p) l l r (p); note that Fi(k,p) does not have a color 
structure. This equation is none other than the self- 
consistent equation in Ref. [4]. By comparing Eq. (2.60) 
with the self-consistent equation in Ref. [4], we see that 
F?(k,p) ee t a ri(k,p)/Cf is the vertex function intro- 
duced in Ref. [4, 7] whose momenta are hard and ul- 
trasoft. We see that the color structure of i" 1 " is the same 
as that of the bare vertex, t a ji, as was shown in Ref. [4] 



-< — ^nnnp — 



¥ 
O 



b ii ii ^ nun -|- 





FIG. 6: The diagrammatic representation of self-consistent equation for Kf at the leading order. For simplicity, A a ^ is set to 
be zero. The solid (curly) line with a black blob is the resummed quark (gluon) propagator, whose information is reflected in 
Eqs. (2.48) and (2.49). 



in a diagrammatic way, by using the fact that r t does 
not have a color structure. 

Also the induced fcrmion source can be expressed in 
terms of the vertex correction: From the definition of 
ilind(X) in Eq. (2.11) and its Fourier transformation, we 
have 



(2 - 61) 



91 



Ki(k,p). 



(27T) 4 

By using the relation [17-19] 

7/ind ( P ) = z r (pMp) 



(2.62) 



which is valid in the linear response regime, and 
Eq. (2.59), we obtain 



d 4 k 
f 

xri(k, P ). 



P(k)pT(k)(nB(k )+n F (k )) 
2k ■ p + 2i(k° + Sm 2 



(2.63) 



This expression coincides with that in Ref. [4]. Thus we 
see that our kinetic equation is equivalent to the self- 
consistent equation in the resummed perturbation [4]. 
This equivalence establishes the foundation of the re- 
summed perturbation scheme. 

By using the correspondence between the self- 
consistent equation in the resummed perturbation the- 
ory [4, 7] and the GBE (Eq. (2.55)), we can obtain ki- 
netic interpretations of the procedures of the resummed 
perturbation theory. Since the structure of the GBE is 
almost the same as that in the QED case [18], we re- 
frain from giving the interpretation of the resummation 
scheme in this paper. 



E. Ward-Takahashi identity 

We show explicitly that the Ward-Takahashi (WT) 
identity derived from SU(N) gauge symmetry is satis- 
fied in the present analysis in this subsection. We also 
show that the WT identity can be derived from the con- 
servation law of the current. 

By multiplying Eq. (2.50) by k^, we obtain 

k^E R ^(k,X) = gt a Vind (X). (2.64) 

See Appendix B for the detailed calculation. This equa- 
tion generates the WT identity at the leading order by 
setting = 0. 

Now we show that the WT identity can be derived 
from the conservation law of the color current, which 
reads [17, 19] 

dtftoM = igt a (*(x)t a mnd (x)-r] ind (x)t a V(x)), 

(2.65) 

where we have set = and introduced j^ nd = t a j^ d . 
By differentiating Eq. (2.65) with respect to ^(y), we get 

d^(y, x) = gt a S(x° - y°)5<®(x - y)t a Vind (x). 

(2.66) 

Here we have set x , y° <E C + , and neglected the term 
which is of order g 3 T^P. By multiplying this equation 
by J d 4 s exp(ifc • s) and taking only the leading-order 



terms [18], wc find 

-k tl 3 Rfi {-k,X)=gt a t a r h ^{X). 



(2.67) 



We see that this equation is nothing but Eq. (2.64) by 
using the color structure of E R,i in the linear response 
regime. 



III. ULTRASOFT FERMION MODE 

In this section, we show the existence of the ultrasoft 
mode in QCD, and obtain the expressions for the pole 
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position and the residue of that mode by solving the 
self-consistent equation that determines the off-diagonal 
propagator in the linear response regime, Eq. (2.58). 
We focus on the momentum region p -C g 2 T with 
p = (p° + i£, p), in which we can solve the self-consistent 
equation analytically. 

As a result of the analysis using Eqs. (2.58) and (2.61), 
we find that the quark retarded self-energy with the mo- 
mentum that satisfies p <C g 2 T has the following expres- 
sion: 



Z R {p) = 



Z ' 



(3.1) 



where Y(p) = p°j° 



-p-7/3 and Z = g 2 C f /(16tt 2 A 2 A 2 ), 
with A = g 2 T 2 Cf/(86m 2 ) and the expression of A is ob- 
tained in Eq. (C5): See Appendix C for a detailed deriva- 
tion. Since Z ~ <j> 2 ; S R (p) is much larger than the inverse 
of the free-part of the quark propagator. Thus, the ul- 
trasoft quark retarded propagator at the leading order 
is 



S R (p) 



27*(p) 

~(- 

2 + IpI/3 + < 



7 



P • 7 



7 + P • 7 
p°-|p|/3 + iC 



(3.2) 



This expression implies that there is a fermionic excita- 
tion in the ultrasoft region, which was suggested origi- 
nally in Ref. [4], whose pole position is 



P=T 



i(. 



(3.3) 



This expression means that that excitation's velocity is 
1/3 and damping rate is C ~ g 2 T\n(l/g). The dispersion 
relation of the mode is plotted in Fig. 7 with those of the 
two soft fermionic excitations [6], for comparison. The 
residue of that pole is 



Z 



l6ir 2 C 



n 2\ * 



f 
N 



g 2 T 2 

5 



8ir 2 (N 2 - 1) \6 



N 



1 

2N 



(3.4) 



If we set TV = 3 and Nf = 3, the expression becomes 



z = g' 



49 
48^2 



(3.5) 



We emphasize that these expressions are obtained in this 
paper for the first time in QCD. These results are summa- 
rized in Table I. We note that the expressions are qual- 
itatively the same as those in the Yukawa model and 
QED [7]. 



IV. ?i-POINT FUNCTION (n > 3) 

In this section, we obtain the equation determining the 
n-point function whose external lines consist of a pair 



TABLE I: The dispersion relation, the damping rate, and the 
residue of the ultrasoft fermion mode in QCD. 



dispersion relation 


T|p|/3 


damping rate 


c f + a ~ s'Tinff- 1 


residue 
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FIG. 7: The dispersion relation in the fermionic sector. The 
vertical axis is the energy p°, while the horizontal axis is the 
momentum |p|. The solid (blue) lines corresponds to the nor- 
mal fermion and the antiplasmino while the bold solid (red) 
one to the ultrasoft mode. The dotted lines denote the light 
cone. Since our analysis on the ultrasoft mode is valid only 
for |p| <C g 2 T, the plot for |p| > g 2 T may not have a physical 
meaning. The residue of the antiplasmino becomes exponen- 
tially small for |p| S> gT, so the plot of the antiplasmino does 
not represent physical excitation for |p| 2> gT, either. The 
parameters are set as follows: g = 0.2 and N = 3. 



of quarks and (n — 2) gluons with ultrasoft momenta, 
where n > 3. We also make an order estimate of the 
n-point function by using that equation. We note that 
the equation was obtained by going beyond the linear 
response theory, and it is a merit of the Kadanoff-Baym 
formalism. 

The fermion induced source also generates n-point 
function not only two point function of the quark, which 
was analyzed in Sec. II C. Here we work in the n = 3 case. 
In this case, we use the following relation [9, 17, 18]: 



S^(p-q-r)gSr^p,-q, 
8 2 rjind{p) 



5V{q)8A a ^(r) 

d 4 fc 5 2 K. l (k 1 p) 



37 



2 , 



(4.1) 



(2tt) 4 5&{q)5A a »(r) 



A,4>=0 



Here gSr^(p, —q, —r) is the quark-gluon vertex correc- 
tion function whose momenta (p, —q, and — r) are ultra- 
soft. To determine 5r®(p, —q, —r), we expand the off- 
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diagonal propagator in the following way: 

Ki(k, X) = Ki(k, X) A=0 + 5Ki(k, X) + 0{A 2 W), (4.2) 

where A^(fc, X) A =o is in the linear order in W while 
5Ki(k,X) is so in <P and A^. By collecting the terms 
that are in the linear order in <P and Eq. (2.54) be- 
comes 

(2ik ■ d x + 2iCk° + 8m 2 ) SK t (k,X) 

-2gk-A a (X)t a K t (k,X) A= o 
= g 2 C f ^p°(k)P^(k)(n B (k ) + n F (k )) (4.3) 



d 4 k' k k "/ j +-f k k' j 
(2tt) 4 k~k' 



SK k (k',X) 



Since Ki(k, X)a=o is determined from Eq. (2.58), we 
can determine SKi(k,X) from this equation. By us- 
ing SKi(k,X), 6r®(p,—q,—r) can be obtained from 
Eq. (4.1). The similar analysis can be done also in the 
case of n > 4. 

Here we make an order estimate of the three point 
function. Since k ~ T and dx ~ g 2 T, we see that 
K A=0 ~ g- Y T-^ by using Eq. (2.58). From this es- 
timate and Eq. (4.3), we find that 8K l ~ g- 2 T~H A» . 
Therefore, the vertex correction whose momenta are ul- 
trasoft is estimated as gST^ ~ g~ x . We note that this 
quantity is much larger than the bare vertex, which is of 
order g. Similarly, in the case that n > 4, the n-point 
function is of order g 2 ~ n . 



V. SUMMARY AND CONCLUDING REMARKS 

We have derived the generalized Boltzmann equa- 
tion from the Kadanoff-Baym equation for quark exci- 
tation with an ultrasoft momentum (~ g 2 T) near equi- 
librium, and showed that the equation is equivalent to 
the self-consistent equation obtained in the resummation 
scheme [4] that is used in the analysis of the quark propa- 
gator with an ultrasoft momentum, in the linear response 
regime in QCD at extremely high T. We have obtained 
the expression of the dispersion relation, the damping 
rate, and the strength of the ultrasoft fermion mode for 
the first time in QCD, by solving the generalized and 
linearized Boltzmann equation. We have showed that 
the Ward-Takahashi identity is satisfied in our approach, 
and that identity can be derived from the conservation 
law of the current. We also have obtained the equation 
that determines the n-point function containing a pair of 
quarks and (n — 2) gluon external lines whose momenta 
are ultrasoft. We note that this equation was obtained 
by retaining the nonlinear response caused by the gluon 
average field in the Kadanoff-Baym equation, and the re- 
summed perturbation has not been able to produce the 
equation. 

As we mentioned in Sec. II A, the possible gauge-fixing 



dependence 2 should be checked by adopting gauge-fixing 
which is other than the temporal gauge. Also, though 
the analyses in this paper and the previous one [7] are 
restricted to the Yukawa model, QED, and QCD, it is 
possible to extend the analysis to other fermion-boson 
system. Especially, investigating the fcrmionic spectrum 
with ultrasoft momentum at high temperature in the 
Weinberg-Salam theory is an interesting task since the 
investigation can be relevant to the analysis of the prop- 
erties of the early universe [26]. 
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Appendix A: SMALLNESS OF LONGITUDINAL 
COMPONENT OF OFF-DIAGONAL 
PROPAGATOR 

In this Appendix, we show that the longitudinal com- 
ponent of the off-diagonal propagator, k t Kf(k,X) 1 is 
negligible compared with the transverse component. 

By multiplying Eq. (2.43) by (# + i@ x /2 - S R (k,X)) 
and subtracting Eq. (2.44) from that quantity, we get 

- k l k 3 K"(k, X) = -n Aabll (k, X)K b j{k, X) 

+ g(S<(k, X)t a \k, X) + $fV(k, X)Df a {k, X)). 

(Al) 

Here we have set A®(X) = for simplicity, neglected the 
terms that are of order g 2 T 2 K°-(k, X), and set k — > i 
in Eq. (2.44). From this equation and Eq. (2.45), we see 
that the longitudinal component of the off-diagonal prop- 
agator, k l Kf(k 1 X), has the same order of magnitude as 
g 2 K?(k,X). 



Appendix B: CALCULATION OF 
OFF-DIAGONAL SELF-ENERGY 

We derive Eqs. (2.51) and (2.64) in this Appendix. 
First we derive Eq. (2.51). To this end, we transform 



2 We note that the analysis in the Coulomb gauge and the temporal 
gauge was performed in QED, and we found no differences in the 
ultrasoft fermion propagator in the two gauge-fixing [18]. 
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Eq. (2.50) in the following way: 



E Hai (k,X) ~ 

■9 / * 1 K A k ' X ) 



(2tt) 4 
+ if abc t c g 2 



k ■ k' 
d A k' 



(Bl) 



P? k {k'-k) (fc' -k)i{k' -k) k 



-2k -fc' (fc'O-fc^^k'-kl 2 , 
x \2k' l g l3 - g %3 {k + k') 1 + 2k j g a ]K]{k',X). 

Here we have used k l Kf{k,X) ~ and k 2 ~ fc' 2 ~ 0, 
and neglected the term that is proportional to k l since 
this term does not contribute to Eq. (2.51) due to the 
presence of P?(k). We also have used ~ 0, 

which can be checked by multiplying Eq. (2.45) by fc" from 
the left. Using k 2 ~ fc' 2 ~ again, we get 



S Rai (k,X) ~ 



(27T) 4 

, d 4 fc' 7 fc 



fc ■ fc' " J 

(k' - %(k' - k) k 



CtCg J (2^) 4 2fc ■ fc' ^°' fe (fc'O - fcO) 2 

x [2k' l g lj - g ij {k + k') 1 + 2k 3 g u ]K b (k',X) 

n 2,,a,b -fabc^ [ k3 l l +l 3 k' 1 b , , > 

i) J /(2^) 4 fc-fc' 



j4j,/ fc 

f abc,c2 I a K 7 



(2tt) 4 fc • fc' 



fc fc (|k'| 2 - fc°fc'°) + fc' fe (|k| 2 - fc°fc'°) 



(fc'° - fc ) 2 



K b {k',X). 
(B2) 



In the last line we have used fc 2 ~ fc' 2 ~ and 
k l Kf(k,X) ~ 0, and neglected the term that is propor- 
tional to fc ? as before. We note that in the derivation of 
Eq. (2.51), the term that is proportional to fc" in Eq. (B2) 
can be neglected since it yields higher order contribution 
after being multiplied by ft. By neglecting the term that 
is proportional to ft and using fiK°;{k,X) ~ 0, we see 
that 7 fc fc fc ~ 7°fc° and 7 fc fc' fc ~ 7°fc'° in Eq. (B2). Using 
these properties, it is easily shown that the g 13 term in 
Eq. (B2) vanishes. Thus that equation becomes 



-Rain 2.a,b [ rf4fc/ tif+^k' 1 b 

a ^ X ^9 t t J (27r)4 — Kj (k ,X). 

(B3) 

By substituting this equation into Eq. (2.45), we will get 
Eq. (2.51). 

Next, we derive Eq. (2.64). By multiplying Eq. (2.50) 



by fc^, we get 
k^H R ^(k : X) 

d A k' 
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+ if abc 9 2 



(27T) 



mv-k) , {k'-kw-k) k 



(fc'-fc) 2 (fc ,o -fc°) 2 |k'-k| 2 / 

x [g l3 \2k' ■ k - fc 2 ) - k 3 (k + k') 1 + k l (2k - k') j ] 



x K b (k',X) 



d 4 k' 
d A k' if abc t c j k 



(t a t b - if abc t c )j 3 K b (k' ,X) 

(fc'-fc),(fc'-fc) fe 



(2tt) 4 (fc' - fc) 2 \ Slk (fc'° - fc ) 2 

x [g lj (2k' ■ fc - fc 2 ) + k 3 (k - fc')'] K b {k',X). 



(B4) 



Here we have used k"K l {k', X) ~ and k' l Kf{k', X) ~ 
in the last line. Using fc' 2 ~ and k li Kf{k',X) ~ 
again, we get 

k^~ Ra »(k,X) 

* 9t a Vind(X) - g 2 if abc t c I ^-^K b (k',X) 



+ if abc t c g 2 



d 4 k' 



:1 



(2tt) 4 ' 3 ' 

(fc' - k) k k 3 



(fc'-fc) 2 (fc'°-fc ) 2 
x ((fc'° - fc ) 2 - (fc' - fc) 2 - |k' - k| 2 )] K b (k', X) 
= gt a r, ind {X). 

(B5) 

We note that the contribution from the gluon self- 
interaction partly cancels the contribution from the 
quark-gluon interaction. 

Appendix C: CALCULATION OF RETARDED 
QUARK SELF-ENERGY WITH ULTRASOFT 
MOMENTUM 

We derive Eq. (3.1) in this Appendix. First, we expand 
the off-diagonal propagator as 

Kdk,p) = Iif\k,p) + K i l 1 \k,p). (CI) 

Here Kf\k,p) (if- (fc,p)) is the zeroth (first) order 
term of the expansion in terms of p/(g 2 T). By collect- 
ing the zeroth order terms in Eq. (2.58) after the Fourier 
transformation, we get 

S,n 2 Kf\k,p) = ^"ikjP^k'UnBlk^+n^k )) 



X gCf 1 3 W{p)+g 



(27T) 4 



fc • fc' 



(C2) 
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To solve this self-consistent equation, we assume that the 
off-diagonal propagator at zeroth order has the following 
structure: 



= 9 -^p\k) P T(k)(n B (k»)+nAk»Wnp), 



(C3) 



where A is a constant. By utilizing this assumption, the 
following expression is obtained [7] : 



d 4 k' fcV + l k k' 3 tA q) 



(2tt) 4 k ■ k' 



(C4) 



where we have used that k is on-shell in Eq. (C2) and 
neglected the term that is proportional to ft, which is 
justified since its contribution to Eq. (C2) is negligible 
due to the presence of another |£. By using this expres- 
sion, Eq. (C2) can be solved and the result is 



By utilizing the properties [7] 



Ag »fT npMp), (cs) 



and 



^A^(fc)V^(fc)(n s (fc°) + n F (k )) 

d A k' fcV+y^ 
(2^ k-k> k (fc ' p) 

8 J (27T)* 7 Kk {k > Ph 



we get 



(C9) 



A = 



1 



1 + A 



(C5) 



Collecting the first order terms in Eq. (2.58) yields the 
following equation: 



Wf^-^-'-Wf^m. (C1 o, 



2k ■ pK^°\k,p) + 5m z K\ y >(k,p) 

= g 2 C0p°(k)PT(k)(n B (k o ) + n F (k )) 

d 4 k' T Ai) nJ , 

~ K k ( k 'P>- 



(C6) 



(27T) 



k ■ k' 



Instead of solving this equation, we multiply this equa- 
tion by 7* and integrate over k. Then, the equation be- 
comes 



A±- 4 2k-pfKl°\k,p) + s m *J 



9 2 Cf 



d 4 k 

WT 4 



\k)fpT(k)(n B (k°) + n F (k )) 



K\ \k ,p). 



(27T) 



k-k' 



(C7) 



Thus, we can calculate the induced fermionic source 
using Eq. (2.61): 



^ind(p) 



Y(p)V(p) 
Z 



(Cll) 



Here we have used the fact that the integral of K^°\k,p) 
vanishes because it is odd function of k. By remembering 
Eq. (2.62), we can obtain the quark self-energy with a 
ultrasoft momentum, which is given in Eq. (3.1), from 
the induced fermionic source given in Eq. (Cll). 

We note that the same result can be obtained from the 
rcsummed perturbation theory [4, 7] because that scheme 
is equivalent to the GBE. 
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